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ABSTRACT 

Using the photometric redshifts of galaxies from the Sloan Digital Sky Survey III (SDSS-III), we identify 
132,684 clusters in the redshift range of 0.05 < z < 0.8. Monte Carlo simulations show that the false detection 
rate is less than 6% for the whole sample. The completeness is more than 95% for clusters with a mass of 
M200 > 1.0 X 10^"^ Mq in the redshift range of 0.05 < z < 0.42, while clusters of z > 0.42 are less 
complete and have a biased smaller richness than the real one due to incompleteness of member galaxies. We 
compare our sample with other cluster samples, and find that more than 90% of previously known rich clusters 
of 0.05 < z < 0.42 are matched with clusters in our sample. Richer clusters tend to have more luminous 
brightest cluster galaxies (BCGs). Correlating with X-ray and the Planck data, we show that the cluster richness 
is closely related to the X-ray luminosity, temperature, and Sunyaev-Zel'dovich measurements. Comparison 
of the BCGs with the SDSS luminous red galaxy (LRG) sample shows that 25% of LRGs are BCGs of our 
clusters and 36% of LRGs are cluster member galaxies. In our cluster sample, 63% of BCGs of rpetro < 19.5 
satisfy the SDSS LRG selection criteria. 

Subject headings: galaxies: clusters: general - galaxies: distances and redshifts 




1. INTRODUCTION 

Clusters of galaxies are usually located at the knots of 
the filamentary structures in the universe. They are the 
most massive bound s ystems to trace t he large-scale struc- 
ture (lBahcall| T988: Post man et all 1 9921: ICarlberg et al.ll 19961: 
iBahcall et al.l Tl997). Statistical studies of clusters provide 
very powerful constraint on the cosmological parameters (see 
a review in [Allen et al. 20 ij) by using, e.g., cluster mass 
funct i on jReiprich & Bohringe? '2002'; 'Seliak' '2002'; 'Dahld 
20Q6t iPedersen & D ahle 2007; Rines et al."2007; Wen et al. 

and gas fraction in massive clusters ( Allen et aLj 
Clusters are also important laboratories t o investi- 
gate the evolution of galaxies in dense environment (iDressled 
Il980l: iButcher & Oemled[T978[ [T984I: iGoto et al.l 12001 "^ 
act as natural telesco pe to study lensed high-redshift faint 
background galaxies (Blain et al. 19991; ISmail et al.l 120021: 
iMetcalfe et al. 2003; Santos et al. 20 0j). As a large number of 
clusters have been detected (IKoester et al.ll2007al: IWen et all 
|2009j), the distributions of galaxy clusters are used to detect 
the baryon acoustic oscillation of the universe (lEstrada et al.l 
I2Q09I: iHtitsil 1201 Ql: iHong et al.l l2Q12h . Correlation of back- 
ground objects with a large sample of clusters shows the ef- 
fects of weak lensing an d spectra line absorption (iMyers et al.l 
l2005l:lL^i^r^rZll2^ 

Galaxy clusters have been found from single-band op - 
tical imaging data (e.g.. lAbelll 11958": Postman etal. 1996), 
multicolor photometric data (e.g., Gladders & Yee 2005; 
iKoester et al. 2007 aD an d spectroscop ic redshift surveys (e.g., 
iHuchra & Geller .19821: lYang et al.l I2QQ5). In addition, 
clusters have also been fo und from X-ray surveys (e.g., 
iBohringer etani2QQ0l |2004'). At millimeter wavelength, the 
Sunyaev-Zel'dovich (SZ) effect has recently used to find clus- 
ters, which is insensitive to cluster redshift (Carlstrom et al.l 
I2OOQ) . Hundreds of SZ clusters have been identified 
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(iMarriage et al.ll201 IblPlanck Collab oration et al."201 Ih. 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) of- 
fers an opportunity to produce the largest and most complete 
cluster sample. It provides photometry in five broad bands (u, 
g, r, i, and z) covering 14,000 deg^ and the follow-up spec- 
troscopic observations. The pho tometric data reach a limit of 
r = 22.2 (IStoughton et al.ll2002h wi th the star-galaxy s epara- 
tion rehable to a limit of r = 21.5 (Lupton et al] 1200 1"). The 
spectroscopic survey observes galaxies with an extinction- 
corrected Petros ian magnitude of r < 17.77 for the main 
galaxy sample (iStrauss et al.l |2002|) and r < 19.5 for th e 
luminous red galaxy (LRG) sample (lEisenstein et al.ll200Th . 
Galaxy clusters or groups have been found by using the 
SDSS spect roscopic data (e.g., Merchan & Zandivarez 20^51: 
[Berlind et al . 2006) and the photometric data (e.g.. iGoto et al.l 
I2QQ2I: IKoester et aL.2007al: IWen et al.ll2QQ9l: iHao et al.ll2Q10l: 
ISzabo etal.ll2Qllh . 

For large samples of galaxy clusters, the determination of 
optical richness still has a large uncertainty because of diffi- 
culties in the discrimination of cluster member galaxies from 
background galaxies using photometric data. This causes a 
large scatter when the optical richness is used to represent 
the cluster mass for the constraint of cosmological parame- 
ters (Hykoff et al. 2008; Rozo et al. 2009b; Wen et al. 20lJ) 

In this paper, we improve the method of IWen et al.l ()2009|) 
using photometric redshifts to identify a large sample of 
galaxy clusters up to 2; ~ 0.8. In Section 2, we identify clus- 
ters using a new cluster detection algorithm and determine a 
cluster richness that is closely related to cluster mass. In Sec- 
tion 3, we compare our sample with the previous cluster sam- 
ples from the SDSS, and correlate the cluster richness with X- 
ray and SZ measurements. In Section 4, we study the bright- 
est cluster galaxies (BCGs) and cross-identify them with the 
SDSS LRGs. A summary is presented in Section 5. 

Throughout this paper, we assume a ACDM cosmology, 
taking Hq =100 h km s"^ Mpc"\ with h = 0.72, = 0.3 
and Qa = 0.7. 

2. CLUSTERS IDENTIFIED FROM SDSS-III 

IWen et al.l (l2009h identified 39,668 galaxy clusters from 
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Fig. 1. — Difference between spectroscopic redshift Zs and photometric 
redshift Zp for galaxies of the SDSS-III. 



the SDSS DR6 by the discrimination of member galaxies of 
clusters using photon i etric r edshifts (hereafter photo- zs) of 
galaxies. IWen & Hani (|2011|) improved the method and suc- 
cessfully identified the high-redshift clusters from the deep 
fields of the Canada-France-Hawaii Telescope Wide (CFHT) 
survey, the CHFT Deep survey, the Cosmic Evolution Sur- 
vey, and the Spitzer Wide-area InfraRed Extragalactic sur- 
vey. Here, we follow and improve the algorith m to identify 
clusters from SDSS-III (SDSS Data Release 8: lAihara etall 
[2QTTh . We first use the SDSS data to determine scaling re- 
lations for cluster mass and radius which will be used in the 
cluster detection. 

2.1. SDSS Data 

The galaxy data of 14,000 deg^ are downloaded from the 
SDSS-III database. The photo-2:s, the iiT-corrections, and 
the absolute magnitudes are obtained from the table Pho- 
to z in which photo- ^ s are estimated based on the method 
of ICsabai et aD feOO?*). We remove objects with deblending 
problems and saturated objects using the flags^. 

As shown in Figure [TJ the uncertainties of photo- zs are 
-0.025-0.030 in the redshift range z < 0.45. They be- 
come larger at higher redshifts. In the following analysis, 
we assume that the uncertainty of photo-z, az, increases with 
redshift in the form of = cro(l + z) for all galaxies. 
Here, we discard those galaxies with a large photo-z error 
zEtt > 0.08(1 + z), i.e., about Sa^, which suffer bad pho- 
tometry or contamination of stars. This procedure removes 
20% of objects, most of which are faint objects (r > 21) with 
large photometric errors. 

BCGs are the most luminous members of clusters. Gener- 
ally, the BCGs are elliptical galaxies and have smaller photo- z 
error than other member galaxies. Proper selection of BCGs 
can be helpful to identify clusters and estimate the cluster pa- 
rameters. To get right BCGs, we select those galaxies as BCG 
candidates which have a photo-z error zErr < 0.055(1 + z) 
and a galaxy ellipticity in the r band less than 0.7. 



2.2. Scaling relations for cluster mass and radius 

The cluster mass and radius are two fundamental parame- 
ters for clusters. The widely used are r2oo, the radius within 
which the mean density of a cluster is 200 times of the criti- 
cal density of the universe, and M200, the cluster mass within 
^200. 

We use the SDSS data of known clusters to get the scaling 
relations for r2oo and M2oo- Here, we take the clusters whose 

3 (flags & 0x20) = and (flags & 0x80000) = and ((flags & 
0x400000000000) = or psfmagerrr^ <= 0.20) and ((flags & 0x40000) 
= 0). 
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Fig. 2. — Upper: scaling relation between cluster radius r2oo and total 
luminosity within a radius of 1 Mpc in units of L*. Lower: scaling relation 
between cluster mass M200 and richness Rl*- The solid lines show the best- 
fit relations. Data are derived for clusters compiled by Wen et al. (2010) . 



masses or radii have been estimated by X-r ay or weak-lensing 
methods as compiled by I Wen et alJ (l2010l) . For each cluster, 
we calculate the total luminosities of cluster member candi- 
dates in the SDSS r-band within a radius of 1 Mpc from its 
BCG by summing luminosities of member galaxy candidates 
brighter than M;(z) < -20.5 but fainter than the BCG within 
a photo- z gap of 2: ±0.04(1 + 2:), with a local background sub- 
traction. Here, is the evolution-corrected absolute magni- 
tude in the r band, M!^{z) = Mr(z) -\- Qz, where we adopt 
a passive evolution of Q = 1.62 ( Blanton et al.l l2003). To es- 
timate local b a ckgro und, we follow the method analogous to 
iPopesso et al.l (|2004|) . For each cluster, we divide the annu- 
als between 2 and 4 Mpc from its BCG into 48 sections with 
equal area. Within the same magnitude and photo-z range, 
we calculate the total luminosity in each sector and estimate 
the mean value and its root mean square. The regions with 
luminosities larger than 3 a are discarded and the mean is re- 
calculated to be the local background. We measure the total 
luminosity of cluster member candidates in units of L*, here 
L* is the evolved characteristic luminosity of galaxi es in the 
r ban d, defined as = = 0)10^-^^^ (Blanton et al.1 

I200I . 

In Figure [21 we show the correlation between cluster radius 
r2oo and total luminosity within a radius of 1 Mpc, LiMpc- 
The best fit gives 

logr2oo = (-0.57±0.09) + (0.44 ±0.05) log LiMpc, (1) 
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where r2oo is in units of Mpc and LiMpc is in units of L*. 
Similarly, we also get the total r-band luminosity within the 
radius of r2oo, ^200, in units of L*, with a background sub- 
traction. In this paper, we define the cluster richness as 
Rl^ = I/2oo/^*- We find that cluster mass M200 and cluster 
richness are closely correlated as (lower panel of Figure [2]) 

log M200 = (-1.49 ± 0.05) + (1.17 ± 0.03) log Rl., (2) 

where M200 is in units of 10^^ Mq. Using these scaling re- 
lations, we can estimate the cluster radius and mass from ob- 
servables of Liupc and L200 or Rl^. We have shown in our 
previous work that the measurements of the total luminos- 
ity or richness a re rob ust for adopting different widths of the 
photo- z gap ( Wen et al.i 2009). In the following, we extrapo- 
late these relations to lower richness for cluster identification. 
Based on the estimated richness, we identify clusters with a 
richness threshold of Rl^ > 12. 

2.3. Cluster detection algorithm 

In I Wen et aP (l2009h . we identified a cluster if more than 
eight member galaxies of Mr < — 21 (without evolution cor- 
rection and different from M^{z) in this paper) are found 
within a radius of 0.5 Mpc and a photo- z gap of z ± 0.04(1 + 
z). The richness was estimated to be the number of galaxies 
within a radius of 1 Mpc which may be systematically smaller 
than the true radius for rich clusters but larger than that for 
poor clusters. Note also that the radius of 0.5 Mpc for clus- 
ter identification is not the same as the radius of 1 Mpc for 
richness estimate. 

With the tight scaling r elations shown in Sec tion | 2. 21 we 
can im prove the method of IWen et al.1 (l2009h and lWen & Hani 
(|2011) to identify clusters from the SDSS-III. Our new algo- 
rithm includes following steps. 

1. For each galaxy at a given photometric redshift, we 
count the number of luminous member galaxies of M^{z) < 
—20.5 within a radius of 0.5 Mpc and a photo- z gap of 2; ± 

0.04(1 + 2;),7Vo.5Mpc. 

2. To Re t cluster candidates, we apply the friend-of-friend 
algorithm dHuchra & Gellejfl982|) to the luminous galaxies 
using a linking length of 0.5 Mpc in the transverse direction 
and a photo-z difference of 0.06(1 + z). The linked galaxy 
with the maximum A/'o.sMpc is taken as the temporary center 
of a cluster candidate. If two or more galaxies have the same 
maximum number, the brightest one is taken as the temporary 
central cluster galaxy. Then, we obtain a temporary list of 
cluster candidates. 

3. For each cluster candidate at 2:, we assume that the galax- 
ies of {z) < —20.5 within a radius of 1 Mpc from the tem- 
porary central galaxy and the photo- 2: gap of 2: ± 0.04(1 + z) 
are member galaxies. The cluster redshift z^ is then defined 
to be the median value of the photometric redshifts of the rec- 
ognized "members". 

4. We recognize the BCG of a cluster candidate as the 
brightest galaxy from the BCG candidate sample within a ra- 
dius of 0.5 Mpc from the temporary center and the photo- z 
gap of z ± 0.04(1 + 2:). As a new development in this pa- 
per, we use the BCG as the center of a cluster and calculate 
^iMpc from the estimated redshift. Subsequently, we get r2oo 
using Equation ([T]) and measure Rl^ from L200 within the es- 
timated r2oo for each cluster candidate. Note that the local 
background has been subtracted already here. 

5. We define a galaxy cluster if Rl* > 12, which corre- 
sponds to M200 ~ 0.6 X 10^^ Mq using Equation ©. Note 
that a cluster candidate at low richness may be contaminated 
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Fig. 3. — Redshift distribution of the 132,684 identified clusters in the 
SDSS-III together with those of 13,823 maxBCG clusters within 0.1 < z < 
0.3 (Koester et al. 2007a), 39,668 WHL09 clusters within 0.05 < z < 0.6 
(iWen et a l. 2009), 55,424 Gaussian Mixture Brightest Cluster Galaxy (GM- 
BCG) clusters within 0.1< z < 0.55 (Hao et al. 2010) and 69,173 adaptive 
matched-filter (AMF) clusters within 0.045 < z < 0.78 ( Sz abo et al.i2011i) . 
(A color version of this figure is available in the online journal.) 



by a small number of bright field galaxies with their photo- 
ns seriously overestimated. To avoid such contamination, we 
require the number of galaxies within a radius of r2oo and a 
photo-2: gap of z ± 0.04(1 + z), N200 > 8. This assures that 
the cluster sample above the detection criteria has a hi gh co m- 
plet eness and a low false detection rate (see Sections [231 and 
12.61) . After we get the cluster candidate list, we merge pos- 
sible repeatedly identified clusters which may survive in the 
previous procedure by using the friend-of-friend algorithm. 
If two cluster candidates have photo-z difference less than 
0.06(1 + z) and a projection separation less than r2oo, we 
merge the cluster candidates as one cluster, and the poorer 
one is then removed from the list. 

Finally, we visually inspect the color images on the SDSS 
Web site"^ for all cluster candidates. About 5000 (3.6%) clus- 
ter candidates are obvious contaminations from bad photome- 
tries and are removed. By this procedure, we discovered 68 
lensing systems from the color images ( Wen et al. 201 1). 

After above procedures, we find 132,684 clusters in the red- 
shift of 0.05 < z < 0.8 from the SDSS-III data. All clus- 
ters are listed in Table [T] Figure [3] show the redshift dis- 
tribution of the identified clusters in the S DSS-III which i s 
compare d with those of the maxBCG (Koe ster et al.ll2007"al) . 
WHL09 ( Wen et al. 2009), Gaussian Mixture Brightest Clus- 
ter Galaxy ( GMBCG: Hao et al 20 10) and adaptive matched- 
filter (AMF; ISzabo et al.1 1201 J samples. Our new sample 
has a peak at 2; ^ 0.42 and the number of our clusters is 
nearly two t imes a s the previously largest AMF sample by 
ISzabo et al.l ('2011). The clusters of z < 0.42 seem to be 
a complete sample. Because of the magnitude limit of the 
SDSS photometric data, the member galaxies ofz > 0.42 are 
incomplete for the threshold of M^{z) = —20.5, so that clus- 
ters ofz > 0.42 have a biased smaller richness than the real 
one because of missing members. In the following, we only 
consider clusters of 2: < 0.42 when we compare our richness 
with other richness estimates and X-ray measurements. 



^ http://skyserver.sdss3.org/dr8/en/tools/chart/list.asp 
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TABLE 1 

Clusters of Galaxies Identified from the SDSS-III 



Name 


R.A.BCG 


Decl.BCG 




2^3, BCG 


^BCG 


^"200 




A^200 


Other Catalogs 




(deg) 


(deg) 






(Mpc) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


WHL JOOOOOO.6+321233 


0.00236 






— 1.0000 


14.92 


1.72 


70.63 


24 


Abell 


WrIL JUUUUUz.3H-UjI /Io 


U.UU^j / 


< OQQT7 


U.Ioyt) 


— i.UUUU 


lo.zU 




1 T /I Q 

I /.4o 






WHL J000003.3+311354 


0.01377 


31.23175 


0.5428 


-1.0000 


20.17 


0.87 


14.27 


8 




WHL J000003.5+3 14708 


0.01475 


31.78564 


0.0932 


-1.0000 


15.18 


0.94 


16.97 


9 




WHL J000004.7+022826 


0.01945 


2.47386 


0.4179 


-1.0000 


19.32 


0.95 


13.71 


10 




WHL J000004.9-033248 


0.02024 


-3.54679 


0.5968 


-1.0000 


20.67 


1.00 


19.19 


11 




WHL J000005.5+354610 


0.02303 


35.76957 


0.4762 


-1.0000 


19.59 


0.91 


15.58 


9 




WHL J000006.0+152548 


0.02482 


15.42990 


0.1656 


-1.0000 


16.60 


1.13 


23.53 


19 


maxBCG,WHL09,GMBCG 


WHL J000006.3+221220 


0.02643 


22.20558 


0.3985 


-1.0000 


19.36 


0.84 


12.73 


11 




WHL J000006.6+100648 


0.02755 


10.11333 


0.3676 


-1.0000 


19.07 


0.93 


16.73 


13 




WHL J000006.6+3 15235 


0.02762 


31.87626 


0.2134 


-1.0000 


17.11 


1.18 


28.58 


15 




WHL J000006.6+292129 


0.02765 


29.35813 


0.2489 


-1.0000 


18.13 


0.92 


15.11 


16 




WHL J000007.1-092910 


0.02957 


-9.48607 


0.3332 


-1.0000 


19.11 


0.80 


12.81 


13 


AMF 


WHL J000007.6+ 155003 


0.03177 


15.83424 


0.1436 


0.1528 


15.99 


1.17 


34.11 


27 


Abell,maxBCG,WHL09,GMBCG,AMF 


WHL J000007.7+ 185245 


0.03208 


18.87909 


0.4347 


-1.0000 


18.95 


1.03 


20.47 


12 





Note. — Column 1: Cluster name with J2000 coordinates of cluster; Column 2: R.A. (J2000) of cluster and BCG; Column 3: Decl. (J2000) of cluster 
and BCG; Column 4: photometric redshift of cluster; Column 5: spectroscopic redshift of BCG, —1.0000 means not available; Column 6: r-band magnitude 
of BCG; Column 7: r2oo of cluster (Mpc); Column 8: cluster richness; Column 9: number of me mber galaxy candidates within r2oo; Column 10: othe r 
catalogs containing the cluster: Abell ( Abell 1958; Abell et al. 1989); maxBCG (Koes ter et al.l2007an : WHL09 ( Wen et al. 2009); GMBCG dHao et al.l20TQh : 
AMF (.Szaboetal...2011i) . 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 



2.4. Redshift uncertainty 

Using the SDSS spectroscopic data, we can get the spectro- 
scopic redshift of identified clusters and verify the accuracy 
of their photometric redshifts. The spectroscopic redshift of a 
cluster is taken to be that of its BCG. We find that the BCGs 
of 38,1 16 clusters have spectroscopic redshifts. 

Figure |4] shows the distribution of the difference between 
photometric and spectroscopic redshifts, — ^s,bcg- The 
standard deviations roughly indicate the accuracy of redshift 
estimate for clusters. In each panel, we fit the distribution of 
% — ^s,BCG with a Gaussian function. The systematic offset 
1^1 of the fitting is less than 0.004, and the standard deviation 
(T is less than 0.018. 

2.5. Completeness of cluster detection 

We use Monte Carlo simulation to test the completeness of 
our sample. Mock clusters are simulated with assumptions for 
their distributions and then added to the real data of the SDSS. 
The cluster detection algorithm is applied to the combined 
data to find the input mock clusters. The detection rate of 
mock clusters is considered as the indication of completeness 
of our detection procedure. 

First, we generate a population of halos with masses 
^200 > Q- ^ . ^ ^P ^^ ^0 following a mass function of 
iJenkins et al ] (I200T I) in a cosmology with = 0.3 and 
cTs = 0.9 in the redshift range of 0.05 < z < 0.60. Then, 
we relate the cluster mass to cluster richness (i.e., total lumi- 
nosity in units of L*) by 

logM200 = -1.49 + 1.17 log(C X i?L*,in), (3) 

where C is a scaling factor between input richness i^L*, in and 
output richness i^L*, out of our detection procedure. The dis- 
crepancy between i^L*, in and i^L*, out is caused by the un- 
certainty of galaxy photo- Because the Equation Q is ob- 
tained for the estimated richness based on the SDSS data, 
we adjust this factor C to make output richness consistent 
with Equation Q. We first set C = 1 for Equation © to 
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Fig. 5. — Correlation between input richness and output richness for the 
mock clusters. The solid line shows the best-fit relation. 
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Fig. 6. — Detection rate of the mock clusters as a function of redshift (the 
upper panel) and cluster mass, M200 (the lower panel). 



get input richness. Member galaxies are simulated for each 
halo within a radius of r2oo- The luminosity function and 
density profile of member galaxies are taken as described in 
IWen et all (l2QQ9l) . Here, the luminosities of member galaxies 
are evolution-corrected, as mentioned in Section 12.21 We as- 



sume that the uncertainty of photometric redshift of member 
galaxies follows a Gaussian probability function with a stan- 
dard deviation of a^, but varies with redshift in the form of 
cr^ = 0.03(1 + 

We add the mock clusters into the real SDSS data and apply 
our cluster detection algorithm to measure the output richness. 
FigureO shows the correlation between the input richness and 
the output richness. The best fit gives 



(-0.52 ± 0.16) + (0.81 ± 0.01)i?L* 



(4) 



Therefore, we get the factor C = 1/0.81 for Equation ^ so 
that the output richness i^L*, out can be related to cluster mass 
by Equation Q. 

Subsequently, we use Equation ^ to get a new population 
of input richness. Again, we apply our method to detect the 
mock clusters. A cluster is detected if the output richness 
^L*,out > 12 and N200 ^ 8. The detection rate is given to 
be the number of detected mock clusters divided by the total 
number of added mock clusters. Figure [6l shows the detec- 
tion rate as a function of redshift. The detection rate is nearly 
100% for clusters with masses M200 > 2x10^^ Mq (richness 
^L*,out > 34) up to redshift ofz^ 0.5, and ~75% for clus- 
ters of M200 > 0.6 X 10^"^ Mq (richness i?L*,out > 12) up to 
redshift of z ~ 0.42. The detection rate drops at higher red- 
shifts, because the member galaxy selection of our algorithm 
is consistent up to z ~ 0.42, as mentioned in Section [231 In 
Figure[6l we also show the detection rate as a function of clus- 
ter mass for clusters of z < 0.42. More than 95% of clusters 
of M200 > 1 X 10^"^ M© are detected. 

2.6. False detection rate 

The presence of the large-scale structures makes it possi- 
ble to detect false clusters because of projection effect. We 
perform a Monte Carlo simulation with the real SDSS data to 
estimate the false detection rate. First, we discard the mem- 
ber candidates of identified clusters within a radius of r2oo 
and a photo-z gap of z^ 0.04(1 -f - z). Second, we shuffle the 
data following l^n & Hani (1201 Ih . Our new cluster detection 
algorithm is applied to the shuffled data. After above proce- 
dures, we get the "false clusters" which satisfy the criteria of 
Rl* > 12 and A^200 ^ 8. The false detection rate is calcu- 
lated to be the total number of "false clusters" identified from 
the shuffled data divided by the number of clusters identified 
from the original SDSS data. We show that the false detection 
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Fig. 8. — Distribution of redshift differences between the WHL12 clusters 
and the maxBCG, WHL09, GMBCG, and AMF clusters within a projected 
separation of r2oo- 



rate is 6.0% for a cluster richness of 12, but decreases to less 
than 1% for a richness of 23 (see Figure [7]). 

Note that we have removed some obvious contaminations 
from cluster list in the cluster identification (Section 12. 3t . 
Here, most of the contaminations cannot be found when shuf- 
fling the cataloged galaxy data. Hence, removing the contam- 
inations slightly decreases the false detection rate. 

3. COMPARISON WITH PREVIOUS CLUSTER SAMPLES 

We compare our cluster sample in this paper (hereafter 
WHL12) with previous ones, including the classical A bell 
sample from the Palomar Sky Survey (Abell et al.lll989i). th e 
maxBCG (Koester et al. 2007a), W HL09 et al.ll2009h . 
GMBCG (Hao etal...201 0) and AMF fezabo et al.ll20 11') clus- 
ter samples from the SDSS. We also compare our clusters with 
X-ray and SZ cluster samples. 



3.1. Comparison with the Abell clusters 

We get the Abell cluster sample (lAbefl et al.lll989 >) from the 
NASA/IPAC Extragalactic Database. For the Abell clusters 
without redshift measurements, we take their redshifts from 
the SDSS photo-2: data. There are 1844 Abell clusters with 
redshifts z > 0.05 in the sky coverage of SDSS-III, of which 
1688 (92%) clusters are within a projected separation of r2oo 
and redshift difference of \Az\ < 0.05 (about 2 a of the un- 
certainty of the redshift difference) from the WHL12 clusters. 
The non-matched Abell clusters are relatively poor clust ers or 
proje cted to neighbors of the large scale structure (Wen et al.l 
I2009h . We mark the Abell clusters in Column 10 of Tableffl 

3.2. Comparison with the maxBCG clusters 

iKoester et al.l (l2007allb h developed a "red- sequence cluster 
finder", maxBCG, to detect clusters dominated by red galax- 
ies. They identify clusters with a richness of A^gai > 10, 
where A^gai is the number of galaxies brighter than 0.4L* in 
the i band within r2oo and 2 a of the ridgeline colors. From 
the sky coverage of the SDSS DR4, they obtained a com- 
plete volume-limited sample containing 13,823 clusters in the 
redshift range of 0.1 < z < 0.3. The maxBCG sample 
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Fig. 9. — Upper: matching rate between the WHL12 clusters and the 
maxBCG clusters, in the same sky coverage of the SDSS DR4, as a function 
of redshift. Lower: matching rate between the clusters of 0.1 < z < 0.3 in 
the two samples as a function of cluster richness. The left panels show the 
matching rate of the WHL12 clusters to the maxBCG clusters, and the right 
ones are the matching rate of the maxBCG clusters to the WHL12 clusters. 



is approximately 85% complete for the clusters with masses 
M> 1 X 10^^ Mq. 

We cross-match the WHL12 clusters with the maxBCG 
clusters in the sky coverage of the SDSS DR4. Figure [8] 
shows the redshift difference (A2: = Zp — 2:maxBCG) be- 
tween the WHL12 clusters and the maxBCG clusters (also the 
WHL09, GMBCG, and AMF clusters for discussion below) 
within a projected separation of r2oo. It shows that a value of 
\Az\ < 0.05 is an appropriate selection threshold for cluster 
matching in redshift difference (top left panel of Figure [8]). 
We find that 10,301 (75%) of 13,823 maxBCG clusters are 
matched with the WHL12 clusters within a separation of r2oo 
and a redshift difference of \Az\ < 0.05. Figure [9] shows the 
matching rates between the WHL12 clusters and the maxBCG 
clusters as a function of redshift. We find that 70%-80% 
maxBCG clusters are matched with the WHL12 clusters in 
the redshift range of 0.1 < z < 0.3, while only 30%-40% of 
the WHL12 clusters are matched with the maxBCG clusters 
because our sample is more complete for poor clusters. The 
matching rates increase with cluster richness for both sam- 
ples. As shown in Figure [9l 95% of the maxBCG clusters 
of ^gai > 20 are matched with the WHL12 clusters, and 
85% of the WHL12 clusters of Rl^ > 40 are matched with 
the maxBCG clusters. Those non-matched maxBCG clusters 
mostly have a richness smaller than the criteria of our cluster 
finding algorithm. A small number of the maxBCG clusters 
of A^gai > 20 are missing in our sample. As we have checked, 
some of them are merged in the friend-of-friend algorithm of 
our procedure (steps 2 and 5), and some of them can be found 
with a larger separation. If we match the clusters within a sep- 
aration of 1.5 r2oo, 97% of the maxBCG clusters of TVgai > 20 
can be found in our sample. 

Figure [TOl shows the correlation between the maxBCG rich- 
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Fig. 10. — Comparison of richness between A^gai in the maxBCG and Rl^^, 
in this paper for the matched clusters. 



ness, A/^gai, and the richness of this paper, Rl^ 
matched clusters. 



for the 



3.3. Comparison with the WHL09 clusters 

From the SDSS DR6, IWen et alJ (I2QQ9 ) identified 39,668 
clusters of galaxies in the redshift range 0.05 < z < 0.6. 
Monte Carlo simulations show that the sample is complete up 
to redshift z = 0.42 for clusters with a richness R > 16.7. 
The cluster richness is the number of member candidates 
of Mr < —21 within a radius of 1 Mpc and a photo-z of 
±0.04(1 + a fter background s ubtraction. 

Comparing to lWen etaP (l2009h , we modify the richness es- 
timate and cluster identification criteria. A cluster is identified 
in this paper when the richness Rl^ > 12 and the number of 
member candidates A^200 > 8 within r2oo, which is looser 
than the c riteria of A^gai > 8 within a radius of 0.5 Mpc in 
IWen et all ( 2009). The cluster sample in this paper is more 
complete for clusters of low richness. The completeness is 
4 0% for clusters w ith a mass of 10^^ Mq (richness of ~12) 
in I Wen et aP (l2009l) , while it is ~90% in this paper. 

Matching is applied between two cluster samples in sky 
coverage of the SDSS DR6. We find that 25,929 (65%) of 
39,668 WHL09 clusters are matched with the WHL12 clus- 
ters within a separation of r2oo and a redshift difference of 
\Az\ < 0.05. Figure [TT] shows the matching rate as a function 
of redshift between the WHL12 clusters and WHL09 clusters. 
The matching rates of the WHL12 clusters to WHL09 clus- 
ters decrease from 85% at z ^ 0.1 to 55% at z ~ 0.4. The 
matching rates of the WHL09 clusters to the WHL12 clus- 
ters increase from 20% at z ~ 0.1 to 40% at z ~ 0.4. The 
matching rates vary as a function of cluster richness within 
0.05 <z< 0.42, 90% of the WHL09 clusters of > 15 are 
matched with the WHL12 clusters, 86% of the WHL12 clus- 
ters of i^L* > 40 are matched with the WHL09 clusters. If we 
match the clusters within a separation of 1.5 r2oo. 91% of the 
WHL09 clusters of > 15 are matched. The non-matched 
WHL09 clusters are probably merged in the friend-of-friend 
algorithm of our procedure, or are poorer than the criteria of 
cluster finding in this work. 
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Fig. 11. — Same as Figure [9] but between the WHL09 clusters and the 
WHL12 clusters in the sky coverage of the SDSS DR6. 
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Fig. 12.— Same as Figure [TOlbut between R in WHL09 and in this 
paper. 



Figure [12] shows a tight correlation between the WHL09 
richness, R, and the richness of this paper, Rl*, for the 
matched clusters. 

3.4. Comparison with the GMBCG clusters 

iHao et al.l (|2010') presented a GMBCG algorithm, which 
detects clusters by identifying the red sequence and the BCG 
feature. From the SDSS DR7, they found 55,424 clusters in 
the redshift range of 0.1 < z < 0.55. They provide two 
richnesses, the weighted richness which is the total number of 
galaxies brighter than 0.4L* within r2oo weighted by a factor 
from the fitting of color distribution, and the scaled richness 
which is the number of galaxies within r2oo and 2 a of the 
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Fig. 13. — Same as Figure [TT] but between the GMBCG clusters and the 
WHL12 clusters in the sky coverage of the SDSS DRV. The solid lines rep- 
resent the matching rate within a separation of r2oo and a redshift difference 
of \Az\ < O.OS.The dashed lines represent the matching rate within a sepa- 
ration of 1.5r2oo and a redshift difference of \ Az\ < 0.1. 
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Fig. 14. — Same as FigurefTolbut between A^gai in the GMBCG and Rl* 
in this paper. 



ridgeline colors. We adopt the weighted richness as GMBCG 
richness A^gais if available, otherwise the scaled richness. 

We cross-match the WHL12 clusters with the GMBCG 
clusters in the sky coverage of the SDSS DR7. We find that 
25,370 (46%) of 55,424 GMBCG clusters are matched with 
the WHL12 clusters within a separation of r2oo and a redshift 
difference of \Az\ < 0.05 (see Figure [13]). The matching 
rates of the WHL12 clusters to the GMBCG clusters decrease 
from 80% at z 0.1 to 40% Sit z ^ 0.4. The matching rates 
of the GMBCG clusters to the WHL12 clusters are 30%-40% 
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Fig. 15. — Same as Figure [bI but between the AMF clusters and the 
WHL12 clusters in the sky coverage of the SDSS DR6. 



within z < 0.4. Within the redshift range of 0.05 <z< 0.42, 
the matching rates increase with cluster richness for both sam- 
ples. We find that 72% of the GMBCG clusters of A^g^is > 20 
are matched with the WHL12 clusters, 83% of the WHL12 
clusters of > 40 are matched with the GMBCG clusters. 
The fraction of missing rich GMBCG clusters is larger than 
those of the maxBCG and WHL09 clusters. This may be due 
to a broader distribution of redshift difference (see Figure [8]). 
If we match clusters within a separation of 1.5 r2oo and a red- 
shift difference of | Az| < 0.1, 55% of the GMBCG clusters 
are matched with the WHL12 clusters, 85% of the GMBCG 
rich clusters of A^gais > 20 are matched (see dashed lines in 
Figure [T3l). 

We note that a few rich GMBCG clusters are still missing. 
Some of them are merged in the friend-of-friend algorithm of 
our procedure, or have a redshift difference of |A2;| > 0.1 
from the WHL12 clusters. Most of other missing GMBCG 
clusters are poorer than the criteria of cluster finding in this 
work. We note a large scatter between the GMBCG richness 
A^gais and (see FigurefT4l). 

3.5. Comparison with the AMF clusters 

ISzabo et aP (1201 ll) used a modified AMF algorithm to iden- 
tify clusters using photometric redshifts of galaxies. From the 
SDSS DR6, 69,173 clusters are identified in the redshift range 
0.045 < z < 0.78 with a richness of A200 > 20. The richness 
of the AMF clusters is defined as the total luminosity within 
r2oo in units of L*, the same as that in this work. 

We cross-match the WHL12 clusters with the AMF clusters 
in the sky coverage of the SDSS DR6. Within a separation of 
r2oo and a redshift difference of \Az\ < 0.05, 27,966 (40%) 
of 69,173 AMF clusters are matched with the WHL12 clus- 
ters. Figure [15] shows the matching rate as a function of red- 
shift. The matching rates of the WHL12 clusters to the AMF 
clusters decrease from 80% at z ~ 0.1 to 30% at z ~ 0.5. The 
matching rates of the AMF clusters to the WHL12 clusters 
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Fig. 16. — Same as Figure flOl but between A200 in the AMF and i?^.* in 
this paper. 



are constantly ~30%. The matching rates vary with cluster 
richness within 0.05 < z < 0.42, 72% of the AMF clus- 
ters of A200 ^ 40 are matched with the WHL12 clusters, 
74% of the WHL12 clusters of Rl^ > 40 are matched with 
the AMF clusters. We note that many AMF clusters have a 
large redshift difference of \Az\ > 0.05 (see Figure [8]). If 
we match clusters within a separation of 1.5 r2oo and a red- 
shift difference of \Az\ < 0.1, 91% of the AMF rich clus- 
ters of A200 ^ 40 are matched with the WHL12 clusters (see 
dashed lines in Figure [T2]). Those non-matched AMF clusters 
are missing due to friend-of-friend merging of our algorithm, 
large redshift difference from the WHL12 clusters and large 
scatter between two richness estimates. 

Figure [16] shows the correlation between the AMF richness 
A200 and the for the matched clusters. 

3.6. Correlations with X-Ray Measurements 

With luminous member galaxies well discriminated in 
the identification procedure, good correlations between the 
cluster richness with X-ray measurements are expected 
(IWen et al.ll2009h . 

We get X-ray clusters from the Northern ROSAT All-Sky 
(NORAS) Galaxy Cluster Survey (Bohringer et al. 2000) and 
the RQSA T-ESO Flux Limited X- ray (REFLEX) Galaxy clus- 
ter survey (iBohringer et al.ll2004|) . The NORAS sample con- 
tains 378 clusters and the REFLEX contains 447 clusters. The 
ROSAT survey provides X-ray luminosity for detected clus- 
ters. We also get X-ray clusters from the query of the X-ray 
Cluster Database (BAX^), where luminosities for 1028 clus- 
ters and temperatures for 210 clusters are available compiled 
from reference. The ROSAT clusters mostly are massive with 
a luminosity greater than 10^^ erg in the 0.1-2.4 keV 
band, while the BAX sample contains many less massive clus- 
ters with a low luminosity less than 10^^ erg s~^. Previous 
work show s that the luminosi ty in both samples is statistically 
consistent (ISzabo et al.ll201l|) . Here, we combine two sample 
for the correlation with our cluster richness. For the clusters 

^ http://bax.ast.obs-mip.fr/ 
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Fig. 17. — Correlations between the cluster richness Rl^ with the X-ray 
luminosity Lx (the upper panel) and temperature Tx (the bottom panel). 
The solid lines are the best fits for the correlations. The dashed lines in the 
upper panel show the icriog Lx best-fit relation. 

in both ROSAT and BAX samples, we adopt the X-ray lumi- 
nosity from ROSAT data. 

X-ray clusters are matched with the WHL12 clusters within 
a separation of r2oo and a redshift difference of \Az\ < 0.05. 
We find 611 matched X-ray clusters with luminosities mea- 
sured, 84 matched X-ray clusters with temperatures mea- 
sured. Figure [17] shows the correlations between cluster rich- 
ness with X-ray luminosity and temperature for the matched 
clusters. The best fits give 

logLx,44 = (-2.49 ± 0.16) + (1.59 ± 0.09) log i^L*, (5) 
and 

logTx = (-0.43 ± 0.05) + (0.61 ± 0.05) logi^L*, (6) 
where Lx^m refers to X-ray luminosity in the 0.1-2.4 keV 
band in units of 10^^ erg s~^, Tx refers to X-ray temperature 
in units of keV. The tight correlations suggest that the richness 
we estimate is reasonable and statistically reliable. The slopes 
of the Lx-R l* and the T x-Rl* relations are in agreement 
with those of ISzabo et al.l (12011). We measure the scatter in 
X-ray luminosity to the best-fit relation and get a value of 
^\ogLx — 0-38, i.e., cTinLx — 0-87. The scatter is in agree- 
ment with that for the correlation of Lx with the maxBCG 
richness ( Rykoff et al. 2008), but larger than those for the cor- 
relations with the improved maxBCG richness (iRozo et al.l 
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Fig. 18. — Correlations between the cluster richness with the SZ measure- 
ment by Planck. 



l2009bl: iRykoff etaP l201 2h. For a similar comparison, we 
measure the scatter in cl uster mass to the best-fit M200-RL* 
relation given in Section [221 and get criog^soo — 0-21, i.e., 
c^inMsoo = 0.48. The scatter in mass based on our richness 
is also in agreement with that based on the maxBCG richness 
(iRozo et al . 2009a). 

3.7. Correlations with the SZ Measurements 

The SZ effect is the result of high energy electrons interact- 
ing with the cosmic microwave background radiation through 
inverse Compton scattering. It has been detected around rich 
galaxy clusters due to a high temperature of host ionized gas 
dCarlstrom et all [19961 l2000l) . The SZ signal is character- 
ized by the quantity D^^^^^^Ysz = {(Jr/'^eC^) J PdV, where 
Da{z) is the angular distance to a cluster at redshift z, ar is 
the Thomson cross-section, c is the speed of light, rrie is the 
electron rest mass, / PdV = J UekTedV is the integration of 
the pressure of hot ionized gas over the cluster volume. As the 
pressure is related to gravitational p otential, D^^^^Ysz is ex- 
pected to be a proxy of cluster mass ([Bonamente et al.[[2008h . 

Recently, the Planck releases an SZ cluster sample includ- 
ing 189 clusters (Planck Collaboration et al. 2011). We find 
that 71 SZ clusters are located at the sky coverage of the 
SDSS-III. Six of them have redshifts less than 0.05. Other 65 
clusters are matched with the WHL12 clusters, of which 61 
have X-ray measurements. Figure [18] shows the correlations 
between the richness and the SZ measurements by the Planck 
for the 61 clusters, where Isz is the SZ signal measured 
within a radius of 5 rsoo and E{z) = y^I^a + l^m(l + z)^- 
The matched clusters mostly have a richness of Rl* > 50. 
Clearly, there is a positive correlation between cluster rich- 
ness and SZ measurement. 
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Fig. 19. — Upper: evolution of BCG absolute magnitude with cluster red- 
shift in the range of 0.05 < z < 0.42 for six richness bins. The average 
absolute magnitudes and the dispersions are plotted for four redshift ranges. 
The average richnesses are calculated and marked for these redshift ranges. 
The solid lines show the best-fitting. Lower: correlation between BCG ab- 
solute magnitude with cluster richness for six redshift bins. Similarly, the 
average absolute magnitudes and the dispersions are plotted for four richness 
ranges, and the average redshifts are marked for these ranges in the plots. 
(A color version of this figure is available in the online journal.) 



4. BCGS AND LRGS IN THE CLUSTERS 

BCGs are luminous and usually located at the centers of 
clusters. Their properties provide the information of galaxy 
formation in the extreme dense environment. Using our sam- 
ple, we can study the evolution of the BCGs. Moreover, we 
cross-identity the BCGs with the SDSS LRG sample. 



4.1. Magnitude evolution of the BCGs 

The BCGs in rich cluste rs are formed at redshift z > 2 
and evolved passively (e.g. JStott et al.ll2008l). The propertie s 
of BCGs are related to their host clusters (IWen & Hanll201lh . 
Richer clusters tend to have brighter BCGs. Here, we use our 



sample to address the magnitude evolution of the BCGs with 
redshift and richness. 

Figure[T9]shows the evolution of absolute magnitudes of the 
BCGs with cluster redshift in the range of 0.05 < z < 0.42 
for six richness bins and with cluster richness for six redshift 
bins. We get the mean values and fit the correlation (solid 
lines) for each bin. Here, the variations of completeness and 
false detection rate are considered by a weighting factor in 
the fitting. The weight is calculated as the percentage of true 
cluster (i.e., 100%— false detection rate) divided by complete- 
ness for each richness and redshift bin. There are very small 
changes of the average richnesses within the redshift ranges 
(upper panels) or the average redshifts within the richness 
ranges (lower panels), indicating that the changes of BCG 
magnitudes have been well decoupled for redshift and rich- 
ness. The best fits between absolute magnitude and redshift 
give 

Mr = (-22.48 ± 0.01) -(1.74 ± 0.03)z 

for 12 < Rl^ < 15; 
Mr = (-22.63 ± 0.01) -(1.76 ± 0.03)z 

for 15 < i?L* < 20; 
Mr = (-22.76 ± 0.01) -(1.77 ± 0.05)^ 

for 20 < i?L* < 25; 
Mr = (-22.93 ± 0.01) -(1.66 ± 0.05)z 

for 25 < Rl^ < 40; 
Mr = (-23.17 ± 0.03) -(1.58 ± 0.09)z 

for 40 < i?L* < 60; 
Mr = (-23.44 ± 0.04) -(1.35 ± 0.17)z 

for 60 < i?L*. (7) 

Clearly, the BCG are brighter in clusters of higher redshifts. 
As the stellar population in the BCGs was formed at redshift 
z > 2 and becomes old with comic time, the BCGs become 
fainter at lower redshifts. To verify the evolut ion of BCG 
magnitude with redshift, using the AMF clusters (ISzabo et aP 
|2Q 1 ID and the brightest one of three BCG candidates for each 
cluster, we get 

Mr,AMF = (-22.43 ± 0.01) - (1.54 ± 0.04)z, (8) 

which is roughly consistent with the results derived from our 
cluster sample. The BCG magnitud e evolu t ion is in agree- 
ment with that of L* found by Blant on et al.l (l2003h . In addi- 
tion, we find that BCGs are brighter in ric her clusters, consis- 
tent with conclusions in previous works (IWen& Had 120111: 
ISzabo et a l. 2011). From the fitting, the BCGs have an av- 
erage absolute magnitude of —22.48 for cluster richness of 
Rl^ = 12-15 at redshift z = 0, but a brighter magnitude 
of —23.44 for cluster richness of Rl^ > 60. This tendency 
can be clearly shown from the correlation between BCG abso- 
lute magnitude with cluster richness for six redshift bins (the 
lower panel of Figure [191). The evolution of the BCGs mag- 
nitudes with redshift and the dependence of richness can be 
formulated together as, 

Mr = {- 21.25 ± 0.01) - (1.75 ± 0.03)z 

-(1.10 ±0.03) log (9) 

4.2. Cross-identification between BCGs andLRGs 

The SDSS LRGs are selected to produce a volume-limited 
sample of massive galaxies to redshift of 0.5. They are ef- 
ficiently selected by color cuts based on passive evolution 
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Fig. 20. — Color-magnitude and color-color diagrams for BCGs in our 
clusters. The solid lines represent the color cuts of the SDSS LRG selection. 



model at redshifts z > 0.15 (lEisenstein et al.ll200Tl) . Gen- 
erally, LRGs are more massive and located at denser envi- 
ronment than general galaxies. They are most likely BCGs 
or brightest group galaxies. BCGs are luminous and mostly 
have red colors. They ar e likely to be LRGs. H owever, this 
may not always be true. ICrawford et al.l (Il999h showed that 
about 27% BCGs in X-ray bright clusters have optical emis- 
sion lines. They may have blue colors due to re cent star 
form ation (Hicks et al. 2010; Wangetal. 2010; Pipi no et aP 
|2Q11|) . Therefore, many BCGs are not LRGs. 

First, we cross-identify the clusters in our sample with the 
SDSS LRGs to address how many LRGs are member galax- 
ies. We get 112,191 LRGs of ^ > 0.15 from the SDSS spec- 
troscopic data. Cross-matching between LRGs with BCGs in 
our sample, we find that 28,336 (25%) LRGs are BCGs and 
40,039 (36%) LRGs are located within a projected separation 
of r2oo and a redshift difference of 0.05. This suggests that 
most LRGs are located not in our clusters but in an environ- 
ment poorer than our clusters. 

On the other hand, we e stimate the fracti on of the BCGs as 
the SDSS LRGs. Follow lEisenstein et al.l (12001) . we match 
the colors of the BCGs with the LRG selection criteria. Fig- 
ure l20l shows the color-magnitude and color-color diagrams 
for the BCGs in our clusters. The color cy is defined as 

C|| = 0.7{g - r) + 1.2[(r - i) - 0.177]. (10) 

If we only consider the color cuts and do not limit the Pet- 
rosian magnitude rpetro, the BCGs can be classified as LRGs 
for those with the color cy above the line in the left panel of 
Figure [20l and the colors, g — r and r — i, between the lines 
in the right panel. There are 83,740 (66%) of 126,041 BCGs 
in the redshift range z > 0.15 satisfying the color cuts of the 
SDSS LRG selection criteria. When the Petrosian magnitude 
limits of Tpetro < 19.2 for Cut I and rpetro < 19.5 for Cut II 
are considered, we find that 63% of BCGs of rpetro < 19.5 
are LRGs. 

5. SUMMARY 

We identify 132,684 clusters of galaxies in the redshift 
range 0.05 < z < 0.8 using photometric redshifts of galax- 
ies from the SDSS-III. The clusters are recognized for those 
with a richness Rl^ > 12 and a number of member galax- 
ies candidates A^200 > 8 within r2oo. Monte Carlo simula- 
tions show that the false detection rate is less than 6% for the 
whole sample. The cluster detection rate is more than 95% 
for clusters of M200 > 1 x 10^^ Mq in the redshift range of 
0.05 < z < 0.42. We cross-match our cluster sample with 
previous cluster samples. The matching rates are 50-70% 
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for previously known clusters but depend on cluster richness. 
Clusters in our sample match more than 90% of previously 
known rich clusters. Comparing with X-ray and the Planck 
data, we find that the determined cluster richness is related to 
the X-ray luminosity, temperature and SZ measurements. The 
richer clusters have brighter BCGs. The BCGs are brighter in 
higher redshift clusters. Cross-matching the BCGs with the 
SDSS LRGs, we find that 25% LRGs are BCGs of our clus- 
ters, 36% LRGs are member galaxies, and 63% of the BCGs 
of Tpetro < 19.5 in our clusters satisfy the SDSS LRG selec- 
tion criteria. 
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